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Abstract 

Manipulation of metal nanoparticles using atomic force microscope is a promising new 
technique for probing tribological properties at the nanoscale. In spite of some advancements 
in experimental investigations, there is no unambiguous theoretical treatment of processes ac- 
companying the movement of metallic nanoislands adsorbed on a flat surface and additional 
research is required. In this paper, we describe computer experiments based on classical 
molecular dynamics in which the behavior of silver and nickel nanoparticles interacting with 
a graphene sheet and sheared with constant force is studied. Frictional force acting on the 
nanoislands is measured as a function of their size. It is shown that its average value grows ap- 
proximately linearly with contact area, and slopes of linear fits are close to the experimentally 
observable ones. The dependence of the friction force value and of the shape of the measured 
friction curves on the type of metal atom is revealed and its possible reasons originating from 
atomistic background are discussed. 

Introduction 

Due to modern trends towards the intensive miniaturization of various devices with moving con- 
stituents friction and wear at the nanoscale is more and more often dealt with in technological 
applications and in many instances these phenomena are the main obstacle to achieving their 
reliable functioning J™ Although traditional experimental techniques, such as friction force mi- 
croscopy (FFM) and surface force apparatus (SFA) have significantly enhanced the understanding 
of atomistic origins of friction, these methods have some deficiencies. In particular, FFM lacks 
the possibility of direct and independent measurement of the true contact area of the sliding in- 
terface, and both mentioned techniques are not able to probe tribological properties of contacts 
with contact area in the range from about several hundreds to some hundred thousands of nm 2 .^A 
promising new approach capable of solving mentioned difficulties is the study of frictional proper- 
ties of adsorbed nanoparticles by moving them with the tip of an atomic force microscope (AFM)J^ 
A variety of experimental works concerned with the manipulation of nanoislands can be found in 
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literature. 5 "'° However, most of them are targeted only at dislocation of the nanoparticles and do 
not investigate tribological processes. Quantified frictional properties of antimony nanoparticles 
grown on highly oriented pyrolytic graphite (HOPG) and pushed with the tip of an AFM have 
been described recently J 5 l 9 l 10 l These experiments indicate the linear dependence of friction force 
on the contact area and the main feature found is that islands with area less than about 10 4 nm 2 are 
much easier to move than the ones with larger contact areas. Additionally, in the last work 9 the au- 
thors revealed the so-called "frictional duality", where some particles with the areas of more than 
10 4 nm 2 assume a state of frictionless sliding while the others show finite friction. However, there 
is no clear explanation of the observed behavior. In earlier works^^ lower friction of the smaller 
nanoparticles is mainly attributed to their compact amorphous structure and incommensurability 
with the substrate in contrast to larger nanoislands, which are often ramified and would not move 
as rigid entities thus providing new route for energy dissipation. While in Ref. 9 the coexistence of 
the two frictional states is ascribed to the presence of contamination molecules and structure of the 
nanoparticles is not considered as a crucial factor defining the observed behavior. 

Mentioned ambiguity indicates the need of additional, in particular, theoretical investigations. 
Existing analytical or seminumerical models™ 111 can provide some estimates of the experimen- 
tally observed quantities. But they are often based on a large amount of assumptions and thus may 
not be able to yield clear picture of the multiple-factor tasks due to neglecting some of the contribu- 
tions which can play appreciable role. Large-scale computer simulations using classical molecular 
dynamics (MD) is an alternative approach capable of elucidation of the detailed atomistic behavior 
of the system and may help to reveal some subtle effects inaccessible for the analytical techniques. 
In literature, however, MD simulations pertaining to the movement of metallic nanoislands on a 
graphite surface are mainly concerned with their diffusion and thus nanoclusters consisting of up 
to several hundreds of atoms are investigated.^^ Such small systems are not appropriate for the 
study of tribological properties and nanoparticles of several orders of magnitude larger should be 
inspected. 

The absence of computational investigations of friction of relatively large metallic nanoparti- 
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cles on a graphitic surface provides the impetus for the current study. In this work we report MD 
simulations exploring friction of Ag and Ni nanoparticles with up to 30,000 atoms on a graphene 
sheet. Such sizes of the clusters are smaller than in the tribological experiments described above. 
Nevertheless, they are suitable for the study of friction (they do not diffuse at room temperature) 
and understanding of their behavior may gain some valuable insights into tribological processes. 
The choice of the metals is caused by the absence of the reliable interaction potential for antimony. 
Properties of Ag (lattice constant, mass, density) are the closest to the ones of Sb among metals 
with face-centered cubic (FCC) lattice^ available from the database of the used potential form. 
Ni is chosen as a trial material for the exploring the influence of type of metal on the system be- 
havior. Additionally, novel experimental techniques of synthesis of Ag-21 and Nf-^ nanoparticles 
with sizes comparable with our model are rapidly developed making these materials possible can- 
didates for the future manipulation experiments allowing direct comparison with the results of the 
simulations. The use of a single graphene sheet in spite of graphite consisting of multiple graphene 
layers can be considered as the first approximation towards approaching the experimental condi- 
tions. Deeper understanding of friction and wear of graphene interacting with various materials 
may also be valuable, as these processes are often related to graphene production J^^l Moreover, 
in literature there is a growing interest directed at the interactions of a graphene sheet with differ- 
ent nanoobjects as they can change electronic properties and structure of this material, which may 
have an implication for future nanodevices.^^The main objectives of our study are to define the 
influence of the size, structure of a nanoparticle and of the type of metal on the friction force and 
to elucidate atomistic processes occurring during the shear of the nanoislands in the context of our 
model. 

Model 

Ideal vacuum conditions are maintained, so our simulations are not able to reveal the role of con- 
taminations in friction of nanoparticles. We consider a graphene sheet lying in the xy plane with 
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zigzag and armchair edges parallel to x and y directions, respectively (see [figure] [!][]![ all snap- 
shots in this work were produced with Visual Molecular Dynamics software' 2 ^'). To hold the sample 
in space, boundary carbon atoms along the perimeter of the graphene layer are held fixed through- 
out the simulations. Silver and nickel nanoislands containing from 5000 up to 30,000 atoms are 
considered. For each nanoparticle's size the unique xxy dimensions of the graphene sheet are used 
and they vary from about 19.68 nmx 17.04 nm to 36.40 nmx31.52 nm, respectively. The total 
number of atoms involved in the calculations varies from 17,800 to 73,808. 

Interactions between carbon atoms in graphene are described by the harmonic potential.^ The 
simulations are not concerned with the formation or breaking of chemical bonds in the layer, so 
this potential form may be more appropriate for our task in comparison with more sophisticated 
potentials such as Brenner 18 19 26 or ReaxFF^I as it is less time-consuming. Forces between metal 
atoms are derived from the alloy form of the embedded atom method (EAM) potential^ which 
is well fitted to basic material properties. For the metal-carbon interaction the pairwise 6-12 
Lennard-Jones (LJ) potential is employed. The values of the parameters e = 0.8738 ■ 10 -2 eV, o = 
2.4945 A with the cutoff distance r c = 2.5c = 6.2363 A from the world^are used for both metals. 
Such an inaccurate approach to the choice of values of the parameters is caused by the absence 
of the reliable data for the interaction energies of Ag and Ni atoms with graphene in literature. 
Although there exist a number of results obtained from the first-principles calculations ^ 1 ' 28 1 29 1 this 
data depends on the relative alignment of the atoms. But in the present problem since metal-C pairs 
conform in many different ways the use of mentioned data would not give the exact interaction 
energy. Moreover, ab initio techniques can overestimate values in several times .^Therefore, our 
approach is a good starting point which should provide a qualitatively correct but not quantitatively 
precise description of the system which should be relevant for metal nanoparticles that weekly bind 
with graphene.^ 

The simulation code is implemented using NVIDIA® CUD A™ platform^HU which allowed 
us to carry out the computations on a single graphics processing unit (GPU) NVIDIA® GeForce™ 
GTX 260. Algorithms for GPU based on the neighbor-list technique from Ref. 31 with our own 
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algorithm for binning atoms into cells are employed, the equations of motion are integrated using 
the leapfrog method^^lwith a time step At = 0.2 fs. 



Results and discussion 

In the course of the simulations such quantities as temperature T of the system, lateral position Xqm 
and velocity of the center of mass (CM) Vx of a nanoparticle, total shear force Fs are measured. 
Also frictional force Ff is defined as the sum of x-components of forces acting on metal atoms 



from carbon ones. Dimensions Lx (cf. [figure] [1][]1 ), Ly, Lz of a nanoparticle are computed as 



the difference between the coordinates of metal atoms with maximum and minimum values along 
the corresponding direction. Structure of the nanoislands is characterized by radial distribution 
function^ (RDF) g(r) computed for time intervals of 1 fs. Although RDF is an integral charac- 
teristic and is not able to reveal the local ordering and subtle structural transitions in the medium, 
nevertheless it is valuable for general characteristic of structure of materials. 

In the experiments^^ nanoparticles are obtained by vapor deposition. Our code allowed 
us to deposit metal atoms on a graphene sheet, but this approach is extremely time-consuming 
for the most of our system sizes. To avoid such a problem in this work the nanoparticles are 
obtained by the procedure imitating the dewetting of thin metallic films by thermal treatments^ 
and in more details it is as follows (movies Ni_16000_tot_35136.avi, Ag_25000_tot_63080.avi, 
Ag_29000_tot_72800.avi are available). At the beginning of a simulation a slab consisting of sev- 
eral layers of metal atoms packed in the ideal FCC lattice is placed above the graphene layer. The 
vertical distance between the graphene plane and the lowest metal layer is 2.1 A and 2.4 A for Ni 
and Ag, respectively. Since many atoms in the slab are on its surface and have coordination which 
is smaller than the one in the bulk state, such ideal FCC lattice is not energetically profitable and 
metal atoms begin to rearrange into the more compact conformation corresponding to the minimal 



free energy. In [figure] [2] [] 2 this is manifested in the decrease of the lateral size of the slab Lx 



with time. This process is accompanied by the release of the energy and therefore the temperature 
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T of the system is being raised (see the bottom plot in |[figure] [2][]2[ ) since at this stage the system 
is in microcanonical conditions, and the nanoisland melts. However, due to much smaller metal-C 
interaction energy comparatively with the metal-metal one, the contact angle of the forming metal 
cluster should approach 180°. 21 The configuration with minimal energy will correspond to a ball, 
which is not suitable for our problem as it has small contact area. So, in order to obtain nanopar- 



ticles with the desired semispherical shape, at the appropriate time moment (in [figure] [2] [] 2 this 
is 26 ps) defined empirically for each system size we begin to apply Berendsen thermostat 26 both 
to metal and graphene during suitable time interval to cool the system down to the temperature of 
about 300 K. After that the thermostat is decoupled from the metal atoms and is applied only to 
graphene to dissipate the heat generated during the shear of the nanoparticle. 

Procedure similar to the described above has been reported in Ref.^ where the formation of 
metal nanoclusters on a graphene sheet for smaller system of 2500 metal atoms was investigated 
using classical MD. However, the obtained in Ref.^2 results showed the formation of a bundle of 
small clusters or even more or less continuous layers of atoms instead of one nanoparticle. Two 
main reasons may be argued for this. Firstly, very low constant temperature of 50 K maintained 
during the simulations might not allow dewetting of the metal atoms. Second reason may be the use 
of a pairwise interaction potential (instead of a many-body one) for metal atoms which is generally 
known not to be suitable for the description of metallic systems.^ 

In Refs.^^ it was found that small nanoparticles (with up to several hundreds atoms) are 
very mobile in the temperature range 300-900 K and diffuse easily on a graphite surface. More- 
over, in Ref.^the authors state, that assuming a proper scaling law for the dependence on size of 
the diffusivity, larger clusters containing up to 25,000 atoms should exhibit significant mobility. 
We decided to verify this, because high mobility of larger nanoclusters means that the energy of 
thermal fluctuations of atoms is comparable with adhesive energy and hence implies zero static 
friction. We considered a graphene layer with dimensions 15.74 nmx 13.63 nm containing 8192 
carbon atoms interacting with Ni nanoisland consisting of 5000 atoms. The interaction potentials 
mentioned in previous section are employed, and it is important to note that in our case the mag- 
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nitude of metal-C energy is about 2.5 times smaller than in Ref. 13 The behavior of the system at 
two values of temperature T = 500 K and 700 K is explored for relatively long simulation runs. 
At 500 K the cooling stage spans the complete simulation time of 5 ns (or 25 million time steps). 
For the case of 700 K the cooling stage is 46 ps and the duration of the simulation is 8 ns (40 
million time steps). These runs are not intended to comprehensively investigate the diffusion of 
the nanoparticle, but the aim is to ensure that in the conditions of our model the smallest nanois- 



lands do not diffuse and therefore friction exists. As [figure] [3] [] 3 shows, at T = 500 K the value 
of Xqm has in fact zero value and is not changed significantly during the simulation, indicating 
the absence of diffusion. However, at 700 K the nanoparticle exhibits high mobility manifested in 
rather sharp changes of Xqm with time (movie Ni_5000_diffuse_T700.avi is available). Thus, one 
may conclude that diffusion grows very rapidly with T. Nevertheless, in contrast to Ref.J^our 
results suggest that it is unlikely for clusters containing more than 5000 atoms to be significantly 
mobile at T less than 500 K. At about 300 K maintained during shear in our simulations significant 
adhesion and friction should exist. 

After the cooling stage the shearing force is applied to the formed nanoparticle. In manipulation 
experiments with AFM pushing but not pulling is always involved.'^' To simulate pushing in our 
system shear force is applied along zigzag edge of graphene (which coincides with the x direction) 
to all metal atoms with values of x-coordinates that are smaller than the x-coordinate of CM Xcm 
of the nanoparticle. At first, the force is incremented in steps of 0.0001 pN until the x-component 
of the velocity of CM Vx reaches the value of 3.55 m/s. Then the shear force acting on each atom 
remains constant, and the simulations are held with constant total shear force Fs- In most cases the 
nanoparticle is translated without rotation in the xy plane. 



[figure] [2] [] 2 partly already discussed shows typical time dependencies of the mentioned quan- 



tities for the Ni nanoparticle containing 16,000 atoms. It can be noted that after the beginning of 
shear Vx and Xqm increase linearly and quadratically in time indicating translation with constant 
acceleration corresponding to the constant Fg. Frictional force Ff acting on the nanoparticle, how- 
ever, is not uniform but has a sawtooth shape with growing in time frequency of spikes, which 
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could indicate the stick-slip motion of the nanoparticle. Such a scenario may be the case for mod- 
els where the objects are translated through the spring moved with constant velocity and where 
sawtooth force may lead to the intermittent behavior of XcM-^But in our case XcM.(t) has contin- 
uous parabolic shape due to constant shear force and therefore the nanoparticle moves continuously 
without intermittency. 



To clarify the behavior of Ff let us consider [figure] [4] []4 and [figure] [5] []5 which plot depen- 
dencies of Ff on Xcm for several Ag and Ni nanoparticles. While for the most of Ni nanoislands 
friction force has a sawtooth form and spikes of Ff are more or less regular, for Ag friction force has 
irregular shape. The distance between spikes for Ni fluctuates near the value of the module of the 
lattice vector in graphene equal to 2.46 A 35 which may indicate the influence of the graphene lat- 
tice on the observed behavior. Note that positive values of Ff are observed suggesting that friction 
force changes direction during the nanoisland translation. On the average the amplitude of spikes 



for Ni is larger than for Ag. Time-averaged value of Ff (see [figure] [6] [] 6) grows approximately 
linearly with contact area A. For Ni relatively large data scattering is observed and Ff has lower 
value for Ag nanoparticles. Slopes of linear fits are 2.91 pN/nm 2 and 1.21 pN/nm 2 for Ni and 
Ag, respectively. The latter value is very close to the experimentally obtained 1 .04 pN/nm 2 for Sb 
nanoparticles with areas more than about 10 4 nm 2 sheared in ultrahigh vacuum.^ Note, that in the 
experiments when Sb nanoparticles are translated on different surfaces (graphite and molybdenum 
disulfide) the obtained linear fits have the same slope. In our case for different metals but for the 
same surface the slopes are different. 

The contact area A in the simulations is defined similarly to the experiments using the lat- 
eral dimensions Lx and Ly of the nanoparticle and approximating it with ellipse. Also note, that 
in |[figure] [2][]2| the dimension Lx linearly grows in time when the nanoparticle approaches the 
boundary of the graphene layer (this was observed for all nanoparticles). To avoid the influence 
of this boundary effect on the results, averaging of Ff is performed from the moment when F$ is 
applied up to the beginning of this change of Lx. 



Shear force F$ and shear stress o also increase approximately linearly with A (cf. [figure] [7] []7 ). 



9 



Values of o are on the order of 10 — 10 2 MPa, which is closer to the experimental data obtained for 
large nanoparticles with areas more than about 10 4 ran 2 . Smaller nanoislands in the experiments 
as a rule require shear stresses on the order of 1 MPa and lower.'^ESl 

Analogous sawtooth shape of the lateral force is often observed experimentally for the tip of 



FFM (cf.™^and references in^% In experiments 4 36 positive values of friction force acting on the 
FFM tip with radius of about 15 nm which is compared with the sizes of our nanoparticles were 



also reported. Additionally, although amplitudes of spikes in^ESare somewhat smaller than in our 
results, the average values of Ff are of the same order as obtained in the present simulations. Such 
a shape of Ff is quite good described by the Tomlinson-type models. 5 37 However, there are some 
discrepancies between the models and the experiments. One of the most significant is that these 
models are often based on the assumption that the tip contains a small number of atoms, ranging 
from one to several tens while in the experiments contact interface contains hundreds or thousands 
of atoms. The atoms are in general incoherently aligned with each other, and the periodic behavior 
should have been washed out since one would be averaging over many atoms out of phase with 
each other. The stick-slip motion with atomic periodicity should not be observed experimentally^ 
but it takes place. This phenomena was detected in the first experiment using FFM with tungsten 
tip probing a graphite surface and at first it was suggested that the reason may be in the cleavage 
of a graphite fiake.LUiM] g u t this explanation fell out of favor once researchers started observing 
atomically periodic friction on nonlayered materials where flake formation is impossible, and the 
question is still opened. 

Our simulations show that sawtooth friction force can be expected during nanoparticle manip- 
ulation. However, it should be noted that the considered system differs from the FFM tip in the 
following aspect. It is known that accordingly to "cobblestone model" two factors contribute to 
friction on small scales: proportional to load and adhesive contribution proportional to the contact 
area.^ While in experiments with FFM Ff can be considerably influenced by the external load,^ 
in our system this contribution is negligible. We estimated the gravitational force which could play 
the role of loading force to be about 5.3 • 10 -20 N for the heaviest Ag nanoparticle of 30,000 atoms. 
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Hence this contribution can be safely neglected (as is done in the model) and therefore adhesive 
part completely prevails over the loading one. 

The simulations also indicate that for the same metal-carbon interaction energy the presence 
of the sawtooth friction force and its average value for the given particle size may depend on the 
type of metal. Such a behavior may be caused by the collective atomic effects resulting from 
the combination of structural and elastic properties of the current metal-graphene system. Let us 
analyze possible factors such as structure and dimensions of the nanoparticles and elastic properties 
of graphene which may have impacted the obtained results. 

In many analytical theories and seminumerical models structure of the surfaces is considered 
to be a crucial factor for the friction.^EES Sliding of ordered commensurate surfaces is predicted to 
occur via stick-slip motion^ while for incommensurate disordered surfaces intermittent motion 
should not be observed. Analysis of the structure of the nanoparticles can be carried out using 



RDF measured at different time moments (see [figure] [8] [] 8) and snapshots of the contact surface 



of the nanoparticles in [figure] [9] []9 As can be seen from [figure] [8] [] 8, after the formation of 
the nanoparticle RDF has completely smeared shape for both metals, indicating the disordered 
structure of the nanoislands. After the cooling phase some ordering is observed manifested in 
the formation of the higher first spike corresponding to the nearest neighbor distance in the bulk 
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state of about 2.49 A and 2.89 A for Ni and Ag, respectively, and of additional smaller spikes. 
However, comparatively with ideal bulk crystal small spikes are very smeared which may indicate 
that nanoparticle is amorphous or has a poly crystalline structure. During shear the structure of the 
nanoparticle is not changed significantly (as follows from RDFs not shown here). Visual inspection 



of |[figure] [9][]9| shows that both Ni and Ag contact surfaces are disordered. Thus, we can conclude 
that generally the surfaces of graphene and of the nanoisland for both metals are incommensurate 
and the long-range order of the surfaces may not significantly influence the shape of the friction 
force. However, one can note that the average distance between the nearest neighbors in Ni near 
2.49 A is very close to the value of the lattice vector in graphene equal to 2.46 A and to the distance 
between spikes in |[figure] [4][]4 and [figure] [5] []5 in contrast to Ag that has larger distance. This 



11 



may cause local commensurability and have an impact on the observed behavior. 



In [figure] [9] [] 9 one can see that the bottom part of the nanoparticles is deformed and due to 
wetting effects lateral dimensions do not exactly reflect the contact area. It is a well known fact 
that AFM images are convolutions of the samples and tip geometries, which gives a correct height 
but an overestimated width of the imaged features^^ and this instrumental error should be taken 
into account. Although in the experiments for relatively large Sb nanoparticles'^^ authors state 
that the area obtained using AFM directly corresponds to the contact area and TEM inspection 
of the interface shows flat contact,^ carefulness should be taken when extracting the contact area 
of the smaller nanoparticles from AFM images. Accordingly to the simulations, wetting effects 
due to small metal-surface interaction energy and the deformation of the substrate may cause 
the inconsistency between the measured by contour and the true contact area. This effect plus 
mentioned instrumental error may cause significant deviation of the measured results from the true 
ones. 

Speaking about the dimension effect, as a rule adhesive interactions are assumed to be propor- 
tional to the area of contact.^ However, in our simulations Ni nanoparticles that have smaller con- 
tact area than Ag experience larger friction. This may be attributed to the quantitative differences 
of the atomic structure of the materials leading to different surface energies of the nanoparticles. 
Smaller nearest-neighbor distance in Ni discussed above may cause the larger number of atoms 
to be located on the surface of the nanoparticle comparatively with Ag, and therefore may lead 
to higher surface energy and stronger adhesion. Additionally, metal atoms located upper than the 
bottom layer may give perceptible for such relatively small nanoparticles contribution to the adhe- 
sion as they lie in the range of the action of the LJ potential. Estimates of the interaction energy 
obtained from LJ potential at distance equal to two lattice constants of a metal give the contribu- 
tion of about 6% and 2% of the energy minimum £ for Ni and Ag, respectively. Thus, taking into 
account the disordered structure of the nanoparticle and the strain of graphene, it is probable that 
metal atoms located farther than the surface layer may give sustainable contribution into adhesion. 
And for Ni it should be several times larger than for Ag. Therefore, friction may depend not only 
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on the area but also on the peculiarities of structure of the nanoparticle in the direction normal to 
the surface. This effect is beyond the capabilities of the standard seminumerical models, such as 
Tomlinson or Frenkel-Kontorova^^ where only one atomic layer is considered. 

Lastly, a few words about the contribution of the dynamic behavior of the graphene layer to 
friction. Fixed atoms at the boundaries of the sheet caused the presence of waves in graphene 
reflecting from these atoms^ and they were not completely damped by the thermostat. Although 
these waves exist both for Ag and Ni nanoparticles, smaller mass and structural peculiarities of the 
latter may be the reasons of higher sensitivity to the waves of the force acting on the Ni nanoparti- 
cles comparatively to Ag ones. This question should be elucidated in the future studies. 

Conclusions 

Atomistic approach employed in the present study allowed us to reveal nonuniform, in particular 
the sawtooth shape of the friction force acting on metal nanoparticles. The simulations indicate 
that for the maintained temperature of about 300 K there is no diffusion of the nanoislands with 
the considered dimensions. But the dependence of the diffusion on temperature may be very rapid 
and it should be investigated in future studies. The results also suggest that time-averaged friction 
force grows approximately linearly with area of contact for both metals, although for Ni large 
data scattering exists. The obtained slope of the linear fit for Ag of 1.21 pN/nm 2 is very close 
to the experimental value for large Sb nanoparticles, and calculated values of shear stress also 
correspond to the ones obtained in the experiment for larger nanoparticles. Generally, in spite of 
ultrahigh vacuum conditions, our simulations did not show the regimes with vanishing friction in 
contrast to experiments described in Ref. 9 

It was also shown that the shape and the average value of the friction force may depend on the 
type of material of the nanoisland and these results are discussed in the context of the peculiarities 
of the structure of the nanoclusters and the dynamic behavior of graphene. The adhesion and hence 
the friction of small nanoparticles may depend not only on the contact area but also on the local 
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structure of a material, in particular in the direction normal to the surface. Also our simulations 
suggest that for small nanoparticles the measured by AFM tip shape of the nanoisland might not be 
the reliable source for the definition of the contact area as there could be significant deviation due 
to surface strains and wetting effects. The waves on the surface may influence the friction and this 
question along with others are the subject of further investigations where the impact of the shear 
direction, elasticity of the layer, number of graphene layers, contaminations, interaction energy of 
metal-carbon atoms should be studied. 
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Figure captions 

Figure 1. Snapshots of the formed nanoparticle containing 25,000 Ag atoms: perspective (a) and 
top (b) views. 

Figure 2. Time dependencies of temperature T of the system, lateral position Xqm and velocity 
Vx of the center of mass of the nanoparticle, total shear force F$, frictional force Ff and lateral 
dimension Lx obtained for Ni nanoisland containing 16,000 atoms. 

Figure 3. Time dependencies of lateral position Xqm of the center of mass, obtained for Ni nanois- 
land containing 5000 atoms at 500 K and 700 K without shear. 

Figure 4. Friction force versus the lateral position of the center of mass of the nanoparticles: Ni 
with 5000 and 25,000 atoms (a), Ag with 5000 and 25,000 atoms (b), Ni and Ag with 16,000 atoms 
(c), Ni and Ag with 29,000 atoms (d). Only initial parts of plots are shown for clarity. 

Figure 5. Friction force versus the lateral position of the center of mass of the Ni (a) and Ag (b) 
nanoparticles containing 29,000 atoms. 

Figure 6. Friction force versus contact area calculated for Ni and Ag nanoparticles. 

Figure 7. Shear stress versus contact area calculated for Ni and Ag nanoparticles. Inset: shear 
force versus contact area. 

Figure 8. Radial distribution function obtained at different time moments for the Ni (a) and Ag 
(b) nanoparticles containing 29,000 atoms. Plots for bulk state are obtained using the same EAM 
potential. 

Figure 9. Side (a,c) and bottom (b,d) views of Ni (a,b) and Ag (c,d) nanoparticles containing 
29,000 atoms. 
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Figure 1: Snapshots of the formed nanoparticle containing 25,000 Ag atoms: perspective (a) and 
top (b) views. 
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Figure 2: Time dependencies of temperature T of the system, lateral position Xqm and velocity 
Vx of the center of mass of the nanoparticle, total shear force Fs, frictional force Ff and lateral 
dimension Lx obtained for Ni nanoisland containing 16,000 atoms. 
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Figure 3: Time dependencies of lateral position Xqm of the center of mass, obtained for Ni nanois- 
land containing 5000 atoms at 500 K and 700 K without shear. 
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Figure 4: Friction force versus the lateral position of the center of mass of the nanoparticles: Ni 
with 5000 and 25,000 atoms (a), Ag with 5000 and 25,000 atoms (b), Ni and Ag with 16,000 atoms 
(c), Ni and Ag with 29,000 atoms (d). Only initial parts of plots are shown for clarity. 
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Figure 5: Friction force versus the lateral position of the center of mass of the Ni (a) and Ag (b) 
nanoparticles containing 29,000 atoms. 
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Figure 6: Friction force versus contact area calculated for Ni and Ag nanoparticles. 



21 



160 




40 60 

A (nm 2 ) 

Figure 7: Shear stress versus contact area calculated for Ni and Ag nanoparticles. Inset: shear 
force versus contact area. 
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Figure 8: Radial distribution function obtained at different time moments for the Ni (a) and Ag 
(b) nanoparticles containing 29,000 atoms. Plots for bulk state are obtained using the same EAM 
potential. 
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Figure 9: Side (a,c) and bottom (b,d) views of Ni (a,b) and Ag (c,d) nanoparticles containing 
29,000 atoms. 
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Figure 10: Table of contents image. 
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